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Abstract 

In weak coupling, the spin gap in doped, even, n-leg periodic Hubbard 
ladders reflects the energy to break a pair into separate quasiparticles. Here 
we investigate the structure of the gap within a spin-fluctuation exchange 
approximation. We also calculate the amplitude for removing a singlet pair 
from two lattice sites separated by a distance (£ x ,£ y ), which describes the 
internal structure of a pair. 
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The half-filled two-leg Hubbard ladder has been found to have a spin-gap |]. Fur- 
thermore, near half-filling, when the interchain hopping is less than twice the intrachain 
hopping, holes doped into the system form singlet pairs and a reduced spin gap remains 
[@||. However, as first noted by Sigrist et al. ||], for the t — J two-leg ladder, the spin 
gap in the doped system differs from the spin gap of the half-filled (undoped) system and 
reflects the energy to separate a pair of holes into two quasiparticles. This is most easily 
understood in the strong coupling limit when the rung exchange J± is large compared to 
the interchain exchange coupling and the hopping t. In this case, the ground state of the 
undoped ladder consists of singlet spin states on each rung and the spin gap is associated 
with exciting one of the rung singlets to a triplet. In the doped ground state there are pairs 
of holes occupying various rungs. In this case, besides the triplet excitation of two spins 
on a rung, a singlet pair can be broken into two quasiparticles which are in a triplet spin 
state. The separate derealization reduction of the kinetic energy of the two quasiparticles 
makes it the lowest energy triplet state. Thus the spin gap is set by the sum of the minimum 
quasiparticle energies. Numerical calculations show that this remains true outside the strong 
coupling limit for the two-leg t—J and Hubbard ladders. A possible weak coupling view of 
the spin gap then is that it represents 2|A p | with p = (pf,0) or (p F ,7i). Here p y = or n 
corresponds to the bonding or antibonding states of the 2-leg ladder. If A p had the simple 
d x 2_ y 2 form (A /2)(cosp a; — cosp y ), then |A P | would be non-vanishing for both (p_f,0) and 

(Pf,7t). 

As the number of legs increases, the size of the spin gap in the half-filled insulating case 
decreases. Thus it could turn out that, the lowest triplet excitation of the doped Hubbard or 
t-J models is similar to the undoped system and simply corresponds to a triplet excitation 
of the spin background rather than to breaking a pair. However, for the 4-leg ladder p y = 0, 
±7r/2, 7r and near half-filling, one can excite quasiparticles at (±7r/2, ±7r/2) where the d x 2_ y 2 
gap vanishes. This would suggest that, at least in weak coupling, near half-filling the spin 
gap of the doped 4-leg Hubbard model vanishes even though the insulating 4-leg Hubbard 
system exhibits a spin gap. Proceeding with the same type of argument for the 6-leg ladder 
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where one can probe p y = 0, ±vr/3, ±2n/3, n, the energy associated with creating two 
quasiparticles would be finite and thus there would be a finite spin gap in the doped 6-leg 
ladder set either by the triplet excitation of the spin background or the two quasiparticle 
triplet excitation. 

Here we investigate the structure of the gap in periodic n-leg ladders within a weak 
coupling spin-fluctuation exchange approximation previously used for the two dimensional 
Hubbard model. We are also interested in calculating the internal structure of the pair wave 
function because it can be directly compared with the results of numerical Monte Carlo and 
density matrix renormalization group (DMRG) calculations. 

The Hubbard Hamiltonian is given by 

Here cJ J(7 creates an electron of spin a on the lattice site and = c\j^Cij^. The 
index j = 1, ...,n denotes a given leg of the ladder, while i denotes the position along a leg 
which we took as 32 sites long. In addition, periodic boundary conditions were used. The 
hopping along a leg is t, the hopping between legs is t± and the onsite Coulomb interaction 
is U. We used a spin-fluctuation interaction 
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V(q,icu m ) = -U 2 x(c[,iuj m ) (2) 

with 

X{q,tUJ m ) = - —— ; r 3 

1 - C/xo(q, 

and 

Xo (q, iu>m) = j7 22 7Z TT- (4) 

p lUJ m - ^p+q - E p ) 

Here e p = —2tcosp x — 2t ± cosp y — /i, uj m = 2mirT, N = n x 32 is the number of lattice 
sites, pi is the chemical potential, and / is the usual fermi factor. The effective Coulomb 
interaction U was adjusted to give strong antiferromagnetic spin fluctuations in x(q, iu) m ) 
and the interaction V(q, ioo m ). 



In order to model the relative internal structure of the pairing fluctuations in an n— leg 
ladder we have studied the eigen solutions of the Bethe-Salpeter equation 



Here p = (p, iu n ) and we have chosen T such that the leading eigenvalue is less than unity. 
In the following we have used the bare one-electron Green's functions G~ 1 (p, iuj n ) = iuj n — e p 
in Eq. (|||). While more elaborate conserving calculations are clearly important in estimating 
T c , here we are interested in developing a picture of the internal structure of a pair. We 
want guidance for Monte Carlo and DMRG calculations and insight into the evolution of 
the spin gap behavior of doped n-leg ladders. 

We first consider the 2-leg case. The bonding (p y = 0) and antibonding (p y = n) bands 
e p = — It cosp x ±2tj_ — p are plotted in Fig. 1 for t± = 0.5i and a filling (n) = 1.0. The RPA 
spin susceptibility x(q, 0), given by Eq. ([|) with U = 1.5i is plotted in Fig. 2. The strong 
antiferromagnetic fluctuations at q = (jr, 7r) for a temperature T = O.li are clearly seen in 
Fig. 2. The eigenfunction 0(p, ioj n ) with the biggest eigenvalue at this temperature is plotted 
in Figure 3 as a function of p x for p y = and tt, and u n = ixT [KJ. As expected, 4>(p,i7iT) 
peaks near the fermi surface and has opposite signs on the bonding and antibonding fermi 
points reflecting the d x 2_ y 2-\ike behavior of the bound state of two holes doped into a half- 
filled ladder. As discussed, since only the wave vectors (jO/,0) and (p/,±7r) are allowed, 
the quasiparticle spectrum is gapped. Thus for the 2-leg ladder a finite spin gap remains, 
consistent with the DMRG @ and renormalization group calculations for the doped 2-leg 
Hubbard ladder. 

The amplitude describing the internal structure of a singlet pair is 



Here A(£ x ,£ y ) is the amplitude for removing a singlet pair from two sites separated by 
(£ x ,£ y ). Using the solution of the Bethe-Salpeter equation, we have 



T 



(5) 



A(£ x ,£ y ) = 2 \(c i+exd+iy ^ jl -c itjl c i+lx:j+ly ^)\^). 
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which is shown in Fig. 4. The relative d x 2_ y 2 like structure in which the amplitude at 
(£ x = ±1, £ y = 0) is out of phase with that at (£ x = 0, £ y = ±1) is clearly seen. 

We have carried out similar calculations for n = 4 and 6 leg ladders. Figure 5(a) shows 
the four band p x dispersion for the n = 4 leg ladder and below it in Figure 5(b) the leading 
Bethe-Salpeter eigenfunction 0(p, iu n ) for u n = irT. In this case the eigenfunction has 
peaks with opposite signs at the p y = and p y — ir fermi surface points and vanishes at 
Py = k/2. This latter point corresponds to {p x = ir/2,p y = tc/2) and represents a node 
in the d x 2_ y 2-\ike eigenfunction. The amplitude A(£ x ,£ y ) is plotted in Fig. 6 and one sees 
the characteristic d^^-like pattern with the 45 degree nodes. Similar results for the 6-leg 
ladder are shown in Figures 7 and 8. In this case, allowed values of p y are 0, ±7r/3, ±2tt/3, 
and 7r. Again, large values of 0(p, ircT) occur at the fermi surface for p y = and n. For 
p y = ±7r/3 and ±2n/3, the gap is small but non-vanishing. The amplitude A(£ x ,£ y ) for 
the 6-leg ladder is plotted in Fig. 8. All of these results were for (n) = 1.0, and we believe 
reflect the internal structure of the pair formed when two holes are added to a half-filled 
n-leg ladder. 

Figure 9 shows the dispersion relations and Bethe-Salpeter amplitude 0(p, inT) for a 
2-leg ladder with (n) = 0.85. Here the two peaks in the bonding and antibonding bands 
arise from umklapp processes. Note that the relative sign difference remains. The relative 
internal pair amplitude (£ x ,£ y ) for this case is shown in Fig. 10. 

Within this approach, the Bethe-Salpeter eigenfunction provides a measure of the quasi- 
particle gap and an n-leg ladder allows one to examine this gap for p y = 0, n/n, 2n/n, ir. 
This is illustrated in Figure 11 for n = 6. Here the fermi surface for the n — > oo two- 
dimensional lattice with t±/t = 0.5 is shown as the dashed curve and <fi(p x ,p y , inT) is shown 
versus p x for p y = 0, ±7r/3, ±2n/3 and ±7r. It is clear that within this approximation the 
spin gap of the doped n-leg ladder is a reflection of the d x 2_ y 2-\ike quasiparticle gap. It will 
be interesting to see to what extent Monte Carlo and DMRG results support this picture. 



5 



ACKNOWLEDGMENTS 



The authors gratefully acknowledge support from the National Science Foundation un- 
der Grant No. DMR92-25027. The numerical computations reported in this paper were 
performed at the San Diego Supercomputer Center. 



6 



REFERENCES 

[1] S.R. White, R.M. Noack and D.J. Scalapino, Phys. Rev. Lett. 73, 886, (1994). 

[2] R.M. Noack, S.R. White and D.J. Scalapino, Phys. Rev. Lett. 73, 882, (1994); Euro- 
physics Letters, 30, 163, (1995). 

[3] L. Balents and M.P.A. Fisher, preprint, |cond-mat/ 9503041) . 

[4] M. Sigrist, T.M. Rice and F.C. Zhang, Phys. Rev. B 49, 12058, (1994). 

[5] Here, the normalization of 4>(p,iu n ) is such that jjY,p4> 2 {p) — 1- 



7 



FIGURES 

FIG. 1. Quasiparticle dispersion e p versus p x for p y = and ir. These results are for a two-leg 
ladder with t± = OM and (n) = 1.0. 

FIG. 2. Magnetic susceptibility x(q>0) versus q x for — and 7r. These results are for a 
two-leg ladder with U = 1.5*, i± = 0.5t, (n) = 1.0 and T = O.K. 

FIG. 3. (f)(p, iitT) versus p x for p y = and 7r. These results are for the same parameters as in 
Fig. 2. 

FIG. 4. (a) Pair amplitude A(£ x ,£ y ) versus (£ x ,£ y ) for the same parameters as in Fig. 2. In 
this figure the intersection of the thick solid lines corresponds to the origin of the 2 x 32 lattice, 

(4,^) = (o,o). 

FIG. 5. (a) Quasiparticle dispersion e p , and (b) </>(p,z7rT) versus p x for different values of p y . 
These results are for a four-leg ladder with U = 1.5t, t± = OM, (n) = 1.0, and T = O.K. 

FIG. 6. Pair amplitude A(£ x ,£ y ) versus (£ x ,£ y ) for the same parameters as in Fig. 5. 

FIG. 7. (a) Quasiparticle dispersion e p , and (b) </>(p, mT) versus p x for different values of p y . 
These results are for a six-leg ladder with U = IM, t± = 0M, (n) = 1.0, and T = O.K. 

FIG. 8. Pair amplitude A(£ x ,£ y ) versus (£ x ,£ y ) for the same parameters as in Fig. 7. 

FIG. 9. (a) Quasiparticle dispersion e p , and (b) </>(p,z7rT) versus p x for p y = and ir. These 
results are for a two-leg ladder with U = 2.0t, tj_ = 0M, (n) = 0.85, and T = O.lt. 

FIG. 10. Pair amplitude A(£ x ,£ y ) versus (£ x ,£ y ) for the same parameters as in Fig. 9. 
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FIG. 11. <fi(p, i7rT) (thick solid curves) versus p x for p y = 0, =br/3, ±27r/3, and ±ir. These 
results are for a six-leg ladder with U = 1M, t± = 0M, T = O.lt, and (n) = 1.0. Here, the dashed 
curves indicate the fermi surface of the two dimensional lattice with t± = 0.5t. 
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Fig. 2 
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